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Abstract — In this paper, the quasi-convexity of a sum of 



quadratic fractions in the form J] 



i=l (1 + dix)' 2 



is demonstrated 
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where a and di are strictly positive scalars, when denned on the 
positive real axis R + . It will be shown that this quasi-convexity 
guarantees it has a unique local (and hence global) minimum. 

Indeed, this problem arises when considering the optimization 
of the weighting coefficient in regularized semi-blind channel 
identification problem, and more generally, is of interest in other 
contexts where we combine two different estimation criteria. 

Note that V. Buchoux et.al have noticed by simulations that 
the considered function has no local minima except its unique 
global minimum but this is the first time this result, as well as 
the quasi-convexity of the function is proved theoretically. 

Index Terms — Asymptotic analysis, Channel estimation, Expo- 
nential polynomial, Minimum MSE, Quasi-convexity, Regulariza- 
tion, Semi-blind estimation Asymptotic analysis, Channel estima- 
tion, Exponential polynomial, Minimum MSE, Quasi-convexity, 
Regularization, Semi-blind estimation A 



I. Introduction 

Many parameter estimation techniques use combined crite- 
ria to exploit different features or properties of the considered 
signals and hence improve the estimation performance. Ex- 
amples of such combined techniques include the blind source 
separation (BSS) method in 0], and the blind equalization 
method in [2] where second and higher order statistics based 
criteria are combined to restore the source signals, and the 
channel equalization and offset estimation technique in Q, 
where again two criteria based on two different features of 
the transmit signals are jointly used to improve the receiver 
performance. 

In U, a similar approach is used for channel shorten- 
ing in OFDM systems, and in Q, O, semi -blind channel 
identification methods are considered where data-aided and 
blind techniques are combined together to shorten the training 
sequence while preserving a high channel estimation quality. 

When combining two criteria, one uses a weighting parame- 
ter that needs to be optimized. In [7|, the weighting coefficient 
is optimized in such a way the asymptotic mean square error 
(MSE) of the channel estimate is minimum. 

The latter is shown to be a non-linear function and its 
optimization in |7j is done numerically using a line search 
algorithm. 

In this paper, we demonstrate that the previous asymptotic 
MSE function is quasi-convex which provides a guarantee that 
numerical optimization always leads to the desired optimal 
weighting parameter value. 



Moreover, one can observe that many asymptotic MSE 
functions have similar forms to that in [7| and thus, we believe 
the result given in this paper might be extended and adapted 
to other problems, where an optimal weighting coefficient is 
needed to combine two contrast functions. As an example, we 
can cite the case where a contrast function is linearly combined 
with an MSE criterion. Referring to the work in JS], we can 
easily show that in this case, the expression of the asymptotic 
MSE has the same form as the one described in 0. 

From the mathematics point of view, our work can be 
viewed as a contribution to the study of the roots of real 
exponential polynomials. It should be noted that this issue has 
been studied for general cases in (9), where interesting results 
about the number of roots of real exponential polynomials 
with real frequencies have been presented. Unfortunately, these 
results yield a loose bound for the number of roots of the 
considered exponential polynomial and thus are of no interest 
for our particular case given by (Q~2}. We have provided in 
our work an original proof that takes into consideration the 
specifications of the considered exponential polynomial. 

The paper is organized as follows. Section II summarizes the 
results in J7J and shows that the considered MSE optimization 
problem can be cast into the optimization problem of a sum 
of quadratic fractions of the form: 



^ 1 + CiX 



i=1 (1 + dixY 



Ci > 0, di > 0. 



Section III and IV are completely devoted to the derivation 
and the proof of the quasi-convexity of the asymptotic MSE 
function. In particular, section III contains some basic notions 
and results about quasi-convex functions. Conclusion and final 
remarks are given in section V. 

Notation: Operators H , _1 and Tr denote Hermitian, matrix 
inversion and trace operators. Moreover, the real and imagi- 
nary parts of a complex z are denoted respectively by Re(z) 
and Im(z). 

II. Regularized semi-blind channel estimation 

In many signal processing applications, the major problem 
is to find out how to estimate some parameters at a low cost 
and with a good accuracy. The best estimate we can have is 
obtained by taking into account all the information that we 
can get about the desired parameter. This approach involves 
in general high computational complexity, thus restricting its 
interest to only theoretical issues. 
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Actually, in practice, suboptimal approaches retain only 
one kind of information on which they are based to derive 
a minimization problem with only a single criterion. An 
intermediate approach that is based on linearly combining 
competitive criteria has been recently proposed in many signal 
processing applications. For instance, in [4], the channel short- 
ening is being improved by linearly combining the null tones 
criterion with that of the guard interval. Also in the context 
of estimating sparse parameter vectors, C p (where p < 1) 
quasi-norms are often linearly combined to standard statistical 
criteria, thus allowing to take into account the sparsity of the 
desired solution ifTOl . 

The optimal selection of the regularizing parameter that 
makes the best trade-off between two different criteria is 
important for the considered parameter estimation problem. 
To the best of our knowledge, the minimization of the mean 
square estimation error with respect to the regularizing coeffi- 
cient has been only analysed rigorously in [7], in the context 
of semi-blind channel estimation. 

Since we will heavily rely on the asymptotic expression 
derived in [7|, it may be illuminating to provide a brief 
overview on the regularized semi-blind estimation technique. 

Regularized semi-blind estimation technique combines blind 
and training based criteria. They have been introduced first in 
the context of Single Input Multiple Output (SIMO) systems 
. In this case, if Sk denotes the unit power transmitted signal, 
the vector received by the N receiving antennas y fc is given 
by: 



1=0 



iSk-i + Vfc 



(1) 



where h; is the N x 1 vector of the Z-th tap of the channel 
impulse response and denotes the additive Gaussian noise. 
We assume that each frame is composed of training and data 
period, (see fig [TJ. The training period corresponds to the 



Pilot symbols Data symbols 

Fig. 1. Time-Multiplexed Training Scheme 

transmission of £ known symbols which are often referred to as 
pilots, whereas the data period corresponds to the transmission 
of p data symbols. 

The blind criterion is based on the statistical properties of 
the received signal in the data period and can be put on the 
form: 

min h H Qh (2) 

l|h||=l 

where h is the channel parameter and Q is a matrix that 
depends solely on the statistical properties of the received 
signal. On the other hand, the training based criterion can be 
expressed as; 

min||y-Sh|| 2 (3) 

h 

where y is the received signal and S is a matrix that depends 
on the pilot symbols. In contrast to blind estimation methods, 



training based techniques are more sensitive to noise and entail 
inefficient bandwidth utilization. However, blind methods are 
more complex, estimate the channel only up to a scalar 
ambiguity and are often non-robust to modelization errors (e.g. 
channel order overestimation errors) [fTTj). For these reasons, 
it might be interesting to combine linearly both criteria so 
as to resolve the drawbacks inherent to blind and training 
based techniques. Hence, the semi-blind estimate is the one 
that minimizes: 



min ||y — Sh| 



Aph H Qh 



(4) 



where A > is the regularizing coefficient and p is the 
length of the information sequence. Note that the semi-blind 
approach in © outperforms the blind approach in (|2) and 
the non-blind approach in <[3j only if the regularizing scalar 
A is chosen properly. In particular, this would be the case if 
A is selected in such a way the asymptotic estimation error 
variance is minimized. Our result is also useful to derive a 
relation between the optimal MSE* and the percentage of 
training symbols which can be adjusted to achieve a target 
MSE performance. 

It has been proved in J7J that the trace of the asymptoticQ 
estimation mean-square error (MSE) is proportional to: 



MSE 



oc Tr{(I + A 7 Qr 1 (l + A 2 7 M(h)) (I + A7Q)" 1 } 



where 7 = ?, and A't(h) is a Hermitian matrix that has 
the same row and column space as Q (meaning that if 
Q = UDU H is the eigenvalue decomposition of Q, A4(h) 
writes as jVl(h) = UAU H , where A is a given Hermitian 
matrix.). 

Using the eigenvalue decomposition of Q, it can be easily 
verified that the MSE is proportional to: 



MSE oc 



1 + A 2 7Q ri 
(1 + X^/duY 



(5) 



where eta > (resp. da > 0) denote the diagonal elements of 
A (resp. the non zero diagonal elements of D). 

Note that in [7 1 and [12|, it was noticed by simulations that 
the MSE has a unique local (global) minimum with respect 
to A, but to the best to our knowledge, until now, this result 
has not been proved in any previous work. 

III. QlJASI-CONVEXITY OF THE MEAN SQUARE ERROR 

For the reader convenience, we recall hereafter the definition 
and also some results about quasi-convex functions, (we refer 
the reader to lfT3l for further information). 

Definition 1: A real valued function / is said to be quasi- 
convex if its domain of definition and all its sublevel sets: 

S a = {x e dom/|/(x) < a} 

for a € R, are convex, where dom/ denotes the set over 
which the function / is defined. 

Examples of quasi-convex functions To illustrate this con- 
cept, we provide in fig [2] some examples of quasi-convex 
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Fig. 2. Examples of quasi-convex functions 



functions. As we can see, we note that a concave and also 
a non convex function can be also quasi-convex. 

Like convex functions, quasi-convex functions satisfy a 
modified Jensen inequality which is given by: 

Theorem 1: A function / is quasi-convex if and only if 
dom/ is convex and for all x, y £ dom(/) and < 9 < 1 

/(fe + (l-%)<max{/(x),/(t/)}. 

Clearly, the quasi-convexity generalizes the notion of convex- 
ity in the sense that the class of quasi-convex functions is 
larger than and includes the class of convex functions. Also, in 
most cases, quasi-convex functions inherit the nice properties 
of convex functions including the absence of local minimum 
as stated in the following theorem. 

Theorem 2: Let / be a quasi-convex function. Then every 
local minimum is a global minimum or / is constant in a 
neighborhood of this local minimum. 

Consequently, if a quasi-convex function / is non constant 
over any given interval (which is the case for the sum of 
quadratic functions we consider), then each local minimum 
is also a global minimum. Moreover, this global minimum 
(whenever it exists) is unique for real valued functions. To 
prove the non existence of local minima besides the global 
one, we use often the following second-order condition: 

Theorem 3: Let / be a real function which is twice deriv- 
able. If / satisfies: 

Vc such that /'(c) = 0, f"{c) > 0, 

then, / is quasi-convex, and each local minimum is a global 
minimum. 

Next we state our main result regarding the unimodality of the 
asymptotic MSE then we prove it in the section after. 

Theorem 4: Let Ci, di be two sequences of n £ W strictly 
positive reals. Then the derivative of 



F n (x) 



71 1 , 2 
y~\ 1 + CiX 



-f x (1 + dixY 



(6) 



Asymptotic refers to the case where p 

R — y ■y 



CO, I 



oo and the ratio 



(2) 

has a unique positive zero xq with F„ (xo) > 0. Conse- 
quently, F n (x) is a quasi-convex function when its domain 
of definition is restricted to M. + and hence has a unique local 
(global) minimum on the positive real axis. In the sequel, we 
will omit the index n for notational simplicity so that F n will 
be referred to as F. 

To prove this theorem, we proceed in the following steps. 

• First, we show that the number of positive real values of 
p j s larger or equal than that of F^ k \ where 
denotes the k-th derivative of F. 

• We introduce the function Gk which has the same number 
of zeros as and prove that it converges uniformly 
to Goo, over a compact set that contains all the zeros of 

p(k)_ 

• Then we prove that G^ has a unique positive zero in 
that compact set. 

• By applying Hurwitz theorem lfT4l . we conclude that for 
large values of k, Gk is zero only once and that will be 
also the case of F( k \ 

• Finally, we prove that the second derivative of F is 
strictly positive when evaluated at the zero argument of 
F. Fig [3] illustrates the shape of function F and its first 
and second order derivatives, for n = 3. 




Fig. 3. Function F and its first and second order derivatives 

Next section provides the details of all these steps and their 
proofs. 

IV. Analysis and Properties of F 
A. Closed-form expressions for the derivatives of F 

In this subsection, we provide a closed form expression 
for the fc-th derivative of function F. We also show that the 
number of zeros of the /c-th derivative is increasing with k. 

Lemma 1: The k-th derivative of F(x) (k > 0) can be put 
on the following expression: 

n , _ 



fW(x) - (-1) 



fc+i \ " 



\k+2 



(7) 



f-f (1 + ekx) 

where a^k and b^k are sequences of positive reals given by: 



bi, k = 2k!c l dt 1 



fc-1 

2d; 
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Proof: See Appendix lAl ■ 
Given the previous expressions of F^ k \ we are able to prove 
our first step result concerning the increasing number of zeros 
of F^ k \ We have the following lemma: 

Lemma 2: Let Zk denote the number of zeros of the k-th 
derivative of F given by (01. Then Zk+i > Zk. 

Proof: Let xx,--- ,x z . denote the zeros of the k-th 
derivative F 1 -^ in [0, oo[. Therefore, using Rolle's Theo- 
rem [15], F( fe+1 ' has at least Zk — 1 zeros yi,--- ,Uz k -i 
where Xi < y, < Xi + i,i £ {1, • • • , Z/. — 1}. Since 
lim^-^+oo F^(x) — 0, there exist at least one zero of F 
in [x Zfc ,oo[. Consequently, the number of zeros of F is 
at least equal to Z^ k \ i.e Z k+1 > Z k . ■ 



Function Gk is defined as: 

Gk (x) = (-l) k +iJJ-e k x k+2 F( k \kx) (11) 

One can easily note that over [r m i n ,r max ], Gk{x) has the 
same number of zeros as F^ k > . 

Clearly, the scaling of the variable x by factor k is in- 
troduced to bring back the roots of F^ from the interval 
[fcr min , kr max ] to the finite length interval [T min ,T max ]. The 

multiplicative function in ( fTTT i (i.e \J~^e k x k+2 ) is introduced 
to normalize the coefficient a. L ,k and bi,k and to approximate 
the denominator terms in (0 by exponential functions. 
Substituting F^ by its expression in (0, Gk writes as: 



B. Uniform equivalence of Gk 

In this subsection, we introduce an alternative function Gk 
that has the same number of positive valued zeros as and 
we provide its asymptotic equivalent expression. For that, let 
us start by providing a useful approximation of coefficient aj,& 
that will be used later to build the function Gk- 

The Stirling formulae [16 | provides us an equivalent @ for 

fc!: 

/• \ k 

k\ ~ \f2nk~ ( 

, e 

we can easily show that: 

a t , k ~ V^k'+h-'c^- 1 ( * + I 

V Ci di 



(8) 



We recall that the overall quasi-convexity proof is based 
on studying the zeros of the function F^ k ' as k goes to 
infinity. Actually, one can show^ that these zeros belong to 
the interval [V k in ,V k ay ], where V k in = min and 

. A lower bound for V^ in and an upper 



V k = 



max 



0-i,k 



G{1,— n} 6 



'i.fc 



bound for V^ ax can be easily computed and are given by: 



V k > k-r ■ 

V k < k-r 
"max — "-'max 



(9) 
(10) 



, d% di 1 
where r min = mm — and r max = max h 7T3 - - 

iG{l,---n} ZCi ie{l,---n} Ci Zdi 

The difficulty that we face is that the zeros of FW 

are of 

order k, thus making the analysis of the asymptotic behavior 
of function F^ somehow delicate. To deal with this difficulty, 
another function, which we denote by Gk, and which brings 
back those zeros to a given fixed interval is introduced. 
This function will be studied over the interval of interest 



'mini 'maxj- 



Equivalence here means that 



-> 1 



/2^fc(|) fc - 

3 Outside this interval, all the terms in the sum given by (TTJ have the same 
sign and hence FW cannot be zero. 



G k (x) = \P^e k x k + 2 J2 kKkX ^ 



i (1 + kdix) 

kb 



( kl 



k+2 
X-l 



=i k k + 2 d k+2 M- + 1 

1 kdix 



E 



9i,k( x ) 
\ h itk (x) 



\/ke 



( kbi 
a i,k [ — 



x — 1 ) and 



where = v ^ d , + . fct+2 

h^) = {kh + l) k+2 . 

In the following, we extend the domain of the function Gk 
to the rectangle R of C given by: 



{z = x + iy 1 x e [r min , r max ] , -e < y < e} 



where e is a constant real that will be specified later. Over 
this domain, the asymptotic equivalent of Gk is given by the 
following theorem: From the previously stated lemma, one can 
prove easily the following result: 

Theorem 5: In the rectangle TZ t , Gk converges uniformly 
to Goo given by: 



where V^i = 577. 

Proof: See Appendix [B] ■ 

C. Zeros of the uniform limit of Gk 

In this section, we prove that Goo has a unique positive real 
zero. This is a byproduct of the following theorem: 

Theorem 6: Let a, bi and cti three sequences of n strictly 
positive real scalars. Let / be the function given by: 



(12) 



Then / admits a unique real positive zero. 

Proof: See Appendix ICl ■ 
By defining f(z) = zGoo(^) and applying Theorem [6] we 
conclude that Goo(z) has a unique real positive zero. 
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D. Application of Hurwitz theorem 

To prove that from a certain range of fc G k is zero only 
once at the real positive axis, we will rely on the following 
known result in complex analysis, [14|: 

Theorem 7: Let fk{z) be a sequence of analytic functions 
in a compact C. Assume that f k converges uniformly to / in 
C. Assume also that / has no zeros on the frontier DC of C. 
Then, there exists fco € N such that Vfc > ko, f and f k have 
the same number of zeros in C. 

Applying this theorem, we can deduce that, G k will have a 
unique zero value in lZ e as Goo, where e is chosen so that 
Goo has no zeros on the frontier of lZ e and has no complex 
zeros besides its real positive zero. Since the number of zeros 
of Gfc is increasing with respect to fc, we conclude that all Gk 
and hence all F 1 -^ have only a unique positive zero. 

Let x z be the unique positive zero argument of F^K Since 
is negative in a neighborhood of zero, and F^ has no 
zeros for x < x z , F^ is negative in the interval [0, x z ]. 
Therefore the function F is decreasing in [0, x z ]. 

Since is positive for large value of x, F^ must change 
its sign at x z , and hence it is positive in the interval [a; z ,oo[. 
Consequently, F is increasing in [x z ,oo[. 

To sum up, we have established that in [0, x z ], F is 
decreasing and in [x z ,oo[ F is increasing. This guarantees 
that x z is a minimum for F and hence F^(x z ) > 0. In 
fact, F^ 2 \x z ) is strictly positive, since F^(x z ) = and 
lim F^(x) = means that there exists y z e ]x z ,oo[ such 

x— >oo 

that F {2 \y z ) = and hence F {2 \x z ) ^ (because F^ has 
a unique zero). 

V. Conclusion 

In this paper, we have provided a rigorous proof for the 
quasi-convexity of the asymptotic MSE of the regularized 
semi-blind channel estimate. 

More generally, we have proved that any function given 
by a finite sum of quadratic fractions J^^a ,c,d > is a 
unimodal function over M + . 

For our considered channel estimation problem, the previous 
result guarantees the absence of non-desired local minima 
of the MSE function when optimized with respect to the 
weighting coefficient. 

Appendix A 
Proof of lemmaQ] 

Proof: Lemma Q] can be proved easily by induction on k. 
For k — 1, we have: 



2cjx(l + dix) 2 - 2dj(l + CiX 2 )(l + d,x) 

^ {i + d lX y 

y-v (1 + dix)(2ax(l + dix) — 2di(l + CiX 2 )) 



(1 + d z x) 



1 {l + dixf 



Let k e N* . Assume that the result is true until order k. Hence, 
Fn can be written as: 



(z) = (-i) fc+1 £?l 



+ d iX ) k + 2 



Therefore, 



F (k+1)( s _ f_i\«+l [ V + d i x ) k+2 

{x) ~ { L) (l + d lx ) 2k +4 

(fc + 2)^(1 + djx) k+1 (b hk x - a itk ) 
(1 + d l x) 2k + i 

_ /-i \k+i h.k + (fc + 2)d l a l ^ k - (AH- l)d t bj^x 
~ l ' ^ (l + d % x) k + 3 

2—1 

_ ( i\fc+2 (fc + l)bj, k djx - (bj,fe + (fc + 2)d i a iik ) 
~ [ ' f-i (l + d t x) k + 3 



i=l 



n , 

- (-} ) k+2 \ " ° l - k + lX " a ',fc+l 
~ [ > 4? (l + d lX ) k + 3 

where 6; ife+1 = (fc + l)b ijk di and a;,k+i = (& 4ifc + (fc 
2)dfOi,fc). Since 6,,/. = 2fc!d* ; ~ 1 c i , we get 1 = 2(fc 

Also, 

a»,k+i = h,k + (fc + 2)dia itk 

= 2k\c l d k - 1 + 2k\(k + 2)d k o, (k ~ ' "'' '"' ' 



= 2(fc + ly.ad 



2 ^/d>2 
, / ! i k + 2fk + ld l k-l 



(fc + fc + 1 V 2 c, ( 2d 



= 2(fc + l)Mf fA , fc + 2 * 



2d; 2 c 4 



Appendix B 
Proof of theorem[5] 

In this lemma, we propose to find the uniform-limit function 
for the function G k 

( z ) = y"_ in 

the rectangle lZ f . 

rev ; A^i=l g i k {z) => £ 

For that, we will first begin by finding the uniform limit 
functions of hi jk and g^k- 

Lemma 3: In the rectangle 7Z e , the sequence of functions 
(h>i,k)k converges uniformly to hi >00 given by: 

i 

Also, the sequence of functions (gi, k ) k converges uniformly 
to 3i,oo given by: 

*- = Ci vis + i) {v ^ z l) 



where Voo.i = \im k - 



'°° Oi k d i XX ' 



Proof: 

The uniform convergence of hi >k to ftj )0O is a by-product of 
the following known result: 

Lemma 4: Over a compact set the sequence function (1 + 
-) n converges uniformly to e z . 
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The uniform convergence of g^f. to gi iOC is obtained by using 
the asymptotic equivalent of a,i t k given in dHJ. ■ 
The uniform convergence of h^k to hi >0o and of g^fc to gi -00 
does not ensure the uniform convergence of j^- to . 

Other extra conditions are needed as it will be noticed in 
the following lemma: 

Lemma 5: Let f k and g k denote sequences of continuous 
functions over a compact C. Assume that g k is bounded over 
C away from zero uniformly in fc and in z, i.e there exists a 
constant M such that: 

VfceN,VzeC \g k (z)\ > M. 



Assume also that f k and g k converge uniformly to /oo and 
goo. Then, ^ converges uniformly to ^ over the compact 

Proof: Since /& and g k are continuous, their uniform 
limits /oo and goo are also continuous. Therefore, there exists 
constant reals Mf, M g such that: 



Vz € C, |/oo(2)| < My and | floo (z)| < M g . 

Since for all fc e N, |</fc(z)| > M, we have |<foo(2)| > M 
To prove the uniform convergence of ^ towards it is 



sufficient to prove that sup 

zee 

tends to infinity. We have: 



fk_ 
9k 



converges to zero as k 



sup 

zee 



It 

9k 



— sup 

zee 

< sup 

zee 

1 

1 



fkg- fgk 



9kg 
fkg - fgk 



< 



M 2 



sup \f k g - fg\ + sup \fg - fg k 

zee zee 



< 



M g sup \f k - f\+ M f sup \g k - g\ 

M \ zee zee 

— ► 



which proves that ^ converges uniformly to ■ 
Since \hi^ k (z)\ > 1 over lZ t , h^k satisfies the condition 
of lemma [5] Applying this lemma on the functions gi tk 
and hi t k, we prove that j^- converges uniformly to 



hi 



9i,k( z ) 



Consequently, G k {z) = C^) \ conver g es uniformly over 

K £ to Goo(z) 



, tk,k(z) 
g%,oo{z) 



Appendix C 
Proof of Theorem|6] 

The proof is performed by induction on n. For n = 1, 
the result is straightforward. Let n 6 N* be a given integer, 
and assume that the result holds true for all k < n, and all 
functions / of the form given by ( fTSl i. Assume that there exists 



di, hi and a,; three sequences of n + 1 strictly positive real 
scalars such that the function 

n+l 



f(x) = J2(^-b i )e-^ x 



i=i 

admits more than one positive zero. Let x% be the first smallest 
zero of / on R + , 

Without loss of generality, we can assume that all the on are 
two by two different and that a n +i = min ou. Since f is 

l<i<n+l 

strictly negative in zero and is positive for large values of x, 
f should change its sign at at least one zero. In the following 
we will consider only the case when / changes its sign at x\. 
The other case coule be treated in the same way. Let x-i be 
the second smallest zero of / on R + . Under this condition, 
we distinuish the following cases: 

• / changes its sign at x\ and at X2- 
» f changes its sign only at x\. 
For the both cases, we can prove that the second derivative of 



g m (x) = e( a "+^>f(x) 



= Y,(a t x - b l )e-( a '~ a ^ 1+ ^) x + (a n+1 x - b n+1 ) e 

i=l 

form e N* 

has three zeros. More particularly, we have the following: 

Case 1: / changes its sign at x\ and at X2 

Since /(0) < 0, f(x) < for x < x%. Therefore, for 
x > X2 and in the vicinity of X2, f(x) < for x > x<2,. 
Since f(x) > for x large enough, / should have a third 
zero £3 > X2- 

For all integers to, we note that / and g m have the same 
number of zeros. Using Rolles theorem, it can be proved that 



the derivative of g r , 
by: 



which we denote jil and which is given 



-(oj - a n+ i H )a l x + bi(a>i 

TO 



x e 



-ie 



a n +i H j 

TO 

(13) 

— [a n +ix 

TO 

(14) 
(15) 



has at least three zeros, since gm (x) tends to zero as x tends 
to infinity. 

Also again by using the Rolle's theorem, we conclude that 
the second derivative of g m denoted by gin (x) has at least 
two zeros. 

Case 2: / changes its sign at only one zero In this case, 
we can also prove that the first derivative of g m has three 
zeros. Actually, at X2, the first derivative of g m must be also 
zero, since X2 is a local minimum for / and hence for g m . As 
g m tends to zero when x tends to infinity, has two zeros 
between ]xi,X2[ and ]x2,oo[. Consequently, in total, gm has 
at least three zeros, and therefore, the second derivative of g m 
denoted has at least two zeros. 
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Taking the derivative of ( fT5l >. (a;) writes as: 



1 

9${%) = y2i a i ~ a n+i + — ) 

*■ — ' m 

i=i 

-2aj - bi(ai - a n+ i + 
2 



a 4 (a 4 - a„ + i H )x 

m 



rn' 



-{pti — a n + 1 + ^)x 



_I, 1 , i \ - 

-a n+ ie m H T(a n+ ix - o„+i)e 



Extending the definition domain of jm' to C + = 
{z — x + iy, x > 0}, we note that for every compact in C + , 
gin 1 converges uniformly to g^ given by: 

71 

goo(z) = - On+i) {ai(ai - a n+1 )z - 2a* 

z=l 

~6,(a,-a rl+1 ))e-^- a "+ 1 ^ 



Let C be the contour corresponding to the rectangle 



R e = < x + iy, x € 



inf 

m,i 



2ai + bi(ai 



sup- 

m.i 



ai(ai - a n+1 + ±] 
2a,i + bi(ai - a n+ i) + 

" ,»€ [-e,ej 



a,(Q!i - a n+ i + —) 

e is chosen such that \goo\ is bounded above zero in C and has 
no complex valued zeros. Then referring to Hurwitz theorem, 

(2) 

goo and gfn will have the same number of zeros in R t for 
large enough values of m, which is in contradiction with the 
induction assumption. 
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